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SUMMARY 

With, the assumption that Berthelot*s equation of state accounts 
for molecular size and intermolecular force effects, and that changes 
in the vibrational heat capacities are given by a Planck term, 
expressions are developed for analyzing one-dimensional flows of a 
diatomic gas . 

The special cases of flow through normal and oblique shocks 
in free air at sea level are investigated. It is found that up to 
Mach numbers of 10 the pressure ratio across a normal shock differs 
little from its ideal gas value;, whereas at Mach numbers above 4 
the temperature rise is considerably below and hence the density 
rise is well above that predicted assuming ideal gas behavior. 

The effects of gaseous imperfections on oblique shock flows are 
studied from the standpoint of their influence on the lift and 
pressure drag of a flat plate operating at Mach numbers of 10 and 
20. The influence is found to be 'small. 


INTRODUCTION 

A wide variety of problems in compressible flow has been 
solved on the assumption that air behaves as an ideal diatomic gas. 
This assumption is justified, provided the pressure and temperature 
range of interest is small and near atmospheric. It is an experi- 
mental fact, however, that when air is subjected to large changes 
in state at pressures and temperatures far removed from atmospheric, 
it ceases to obey the simple gas law and exhibits other properties 
not characteristic of an ideal gas. Consequently, flow processes 
in which air is subjected to these extreme conditions can be 
expected to depart from perfect gas behavior. It is known that 
such flowB will be encountered by aircraft flying at high supersonic 
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speeds and in hypersonic wind tunnels; hence, the nature and extent 
of this departure have become important considerations in aerody- 
namics. 

Classical theories and experiments have shown that three 
properties of a real gas first cause it to exhibit characteristics 
unlike those of an ideal gas. These properties may be classified 
as thermal and caloric imperfections. Thermal imperfections in 
the form of intermole cular forces and molecular size effects are 
significantly manifest at low temperatures and high pressures. 

Changes in the vibrational heat capacities become an important 
. caloric imperfection at relatively high temperatures. Circumstances 
under which effects of molecular dissociation and electronic excit- 
ation become important (e.g., temperatures appreciably above 5000° R) 
may be neglected for the present. Insofar as gases in equilibrium 
are concerned, it is usually sufficient to account for intermole c— 

; ular force and molecular size effects with additional terms in the 
equation of state. Similarly, changes in the vibrational heat 
capacities of the molecules may be accounted for with a function 
of temperature in the expressions for the specific heats. 


Tsien (reference' l) investigated the effects of gaseous 
inperfections on air flows using Van der Waals state equation. 
Approximate solutions to the one— dimensional isentropic and normal 
shock equations were obtained. (He points out two very early 
papers of limited scope on the subject by A. Busemann and 
W. J. Walker.) The Joule— Thomson effect was neglected in Tsien's 
analysis 1 , however, thus introducing some Brror. Donaldson 
(reference 2) observed this error and found that the differential 
equations of motion could be integrated to yield one— dimensional 
isentropic flow equations, exact to the accuracy of Van der Waals 
equation. A comparison of results obtained with these equations 
and the flow equations for a perfect gas indicated that appreci- 
able error would be made if the latter expressions were used to 
predict aerodynamic phenomena involving high temperatures or high 
pressures. 

A matter considered by Kantrowitz (reference 3 ) is the 
inability of a gas composed of polyatomic molecules to instan- 
taneously adjust its internal energy to temperature changes at high 

^■This error is reported to have been subsequently corrected; however, 
the author is unaware of the publication. 
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temperatures. This time lag in equilibrium partition of energy 
occurs primarily in the vibrational energies of the molecules 
which, unlike the translational and rotational energies, require 
many collisions per molecule to became fully adjusted to a new 
temperature. As pointed out by Eantrcwitz, for rapidly acceler- 
ating flows the resulting heat— capacity lag causes an entropy 
increase and hence a total-head decrease in a gas. Tha net effect 
may be made negligible, however, by decreasing the rate of change 
of temperature of the fluid with time. In the special case of flow 
through a normal shock wave, Bethe and Teller (reference k) found 
that heat-capacity— lag effects caused the thickness of the shock 
to increase. (The wave is considered to extend over that portion of 
the fluid in which equilibrium partition of energy does not prevail.) 
As the Mach number upstream of the wave increases, however, the 
shock tends rapidly toward a sharp discontinuity. It may also be 
observed that, although heat-capacity lag influences flow within 
the shock, it does not alter the magnitude of the entropy increase 
across the shock , as this is independent of the mechanism by which 
available energy 1 b dissipated within the wave. These observa- 
tions are important for they lend support to the simplifying 
assumptions that moderately accelerating flows may be treated as 
isentropic except in the presence of a shock wave, in which case 
the wave may be treated as a discontinuity and the properties on 
both sides calculated as if there were no heat-capacity lag. 

In the present paper, one-dimensional flow equations are 
obtained with the aid of Berthelot's equation of state. This equa- 
tion, rather than Van der Waals , is employed in order that somewhat 
better estimates of intermolecular force effects may be obtained. 

It Is assumed that calorio imperfections may be accounted for with 
a Planck term in the expressions for the specific heats. The 
analysis covers both isentropic and plane-shock flows. Within the 
limitations of the assumptions, both exact and approximate flow 
equations are developed. The approximate expressions provide 
explicit solutions for all the Important flow parameters, and are 
designed to utilize data available on the theoretical behavior of 
an ideal diatomic gas. (See, e.g., reference 5.) 


SYMBOLS 


a local speed of sound, feet per second 
A nozzle cross-sectional area, square feet 



NACA TN No. l86l 


b molecular size constant, cubic feet per slug 

c intermolecular force constant, °R, feet to the fifth per slug, 

second squared 

c p specific heat at constant pressure,' foot-pounds per slug, °R 
c y specific heat at constant volume, foot-pounds per slug, °R 
e base of natural logarithms, 2.718. . . 

Zn logarithm to base e 

M Mach number (ratio of local velocity to local velocity of 
sound) 

P absolute pressure, pounds per square foot 

R gas constant ^1715 for air^ > foot-pounds per slug °R 

T absolute temperature, °R 

u internal energy, foot-pounds per slug 

V local velocity, feet per second 

v specific volume , cubic feet per slug 

w external work performed, foot-pounds per slug 

7 ratio of specific heats (cp/c v ) 

p mass density, slugs per cubic foot 

8 molecular vibrational energy constant, °R 


Subscripts 


stagnation conditions 


conditions upstream of shock wave 
conditions downstream of shock wave 


i 


ideal gas quantities 
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Superscripts 

* quantities at the nozzle throat 

' quantumr-mechanical functions 

ANALYSIS 

It will be assumed throughout the analysis that the gas exhibits 
no heat-capacity lag. Shock-free flows are considered isentropic, 
and flow through plane shock waves is assumed to be adiabatic. Justif- 
ication for the first two assumptions is discussed in the introduction. 
The extent to which plane shock processes, particularly at high Mach 
numbers, deviate from adiabatic behavior is left for future investiga- 
tion. 


Berthelot’s equation of state is employed in two forms, depend- 
ing upon the desired range of applicability. In the development of 
equations for investigating flows over a wide range of Mach numbers, 
temperatures, and pressures, the state equation is used in its exact 
form, 

p = P ET _ °P 2 

1-bp T (1) 

where b is the molecular size constant and c is the inter- 
molecular force constant. It will be noted that equation (l) differs 

from Van der Waals equation by a factor of i in the intermolecular 

T 

force term. The introduction of this factor yields a variation of 
intermolecular force with temperature which is in accordance with 
experiment. (See, e.g. , reference 6.) Thus, as pointed out in 
reference 7 , close agreement between theoretical and experimental 
values of compressibility of gases is obtained. 

In the derivation of approximate flow equations an expression 
of the form 


P = pRT (1+bp - j^g) (2) 

is employed. This equation is, of course, a first-order approxi- 
mation to equation (l). 


Flow of a Diatomic Gas Obeying Berthelot's Equation of State 

Isentropic flow of a Berthelot gas may be studied with the aid 
of the differential isentropic expansion equation. 
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du + dw = 





dT + Pd 



= 0 


where 



and the Joule-Thorns on effect ia given by 



( 3 ) 


( 4 ) 


Combining these expressions and substituting the value of P from 
equation (l) yields 


RT 


c v dT - j dp - dp = 0 


( 5 ) 


Now the differential expression for c v is 



which, upon substituting from equation (l), may be integrated to give 

c v = c v’ + (6) 

where c v * is a function which describes the quantum-mechanical 
variations of Cy with temperature . The second term on the right 
of equation (6) represents the effects of gaseous imperfections on 

Cy • 


The function chosen for Cy f is determined by the molecular 
structure of the gas under consideration and the temperature range 
over which accurate predictions of c v f are desired. For aerody- 
namic purposes * diatomic gases are of primary interest. The 
important temperature range extends from liquefaction temperatures 
to several thousand degrees Rankine. A relatively simple function 
may be written for c v * in this case, as the number of trans- 
lational and rotational degrees of freedom is constant, and only the 
variation with temperature of the vibrational heat capacity need be 
considered. This function is 
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c 


I 

V 


°Vi 




1— e 


(e/T) 

Wt)12 


(7) 


The second term in the brackets, essentially a Planck term, accounts 
for the vibrational contribution to the specific heat at constant 
volume. 2 The assumption is that the molecules of the gas behave like 
linear harmonic oscillators insofar as the vibrational degrees of 
freedom are concerned. (See reference 8.) 


An expression governing isentropic expansion of an imperfect 
diatomic gas may now be obtained by substituting equations (6) and 
(7) into equation (5) and integrating from stagnation to static 
states. Since 


/p_\ _ T dp — 2pdT 

W ' 


rp3 


there results the relation 


[e (9/T °>-l] 


Cy.ln fifr) +R ln p ^ 1 ~ bp | - c ( -2- - -22-') + E In — 

1 VW p(i-t>p 0 ) VT 2 T 0 V " e ( 0 / T )_ 1 J 


+ B 


0 e (6 / T) 

T 


= 0 


--J T o 


(8) 


In order to determine the Mach number of a stream, it is neces- 
sary to find the velocity of flow and speed of sound in the stream. 
These quantities may be found by employing the one— dimensional 
energy equation. 


du + d(Pv) + VdV = 




dT + d 



+ VdV = 0 


(9) 


This is a common method of accounting for the variation with temper- 
ature of the vibrational heat capacities. It has been adopted by 
Donaldson and others for imperfect gas studies. 
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Substituting equations (l), (4), ( 6 ), and (j) into equation ( 9 ) and 
integrating from stagnation to static temperature and density yields 
for the velocity 

To 

The corresponding speed of sound is determined by substituting equa- 
tions (l), (5), (6), and ( 7 ) into the general equation 



The resulting expression is 


RT _ 2cp 
(1-bp) 2 T 


Combining equations (10) and (ll) yields the following equation for 
the Mach number. 



( 12 ) 


The specific heats and the ratio of specific heats are readily 
obtainable for a diatomic gas obeying Berthelot's state equation. 
The specific heat at constant volume c y is found by substituting 
equation ( 7 ) into equation (6), thus yielding 
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= a 


r ± {l+(n-D 



e (e/T) 

'"l-* (0/T)] 2 + 


2cp 

RT 2 - 


} 


(13) 


The specific heat at constant pressure Cp is obtained by substi- 
tuting equations (l) and ( 13 ) into the reciprocity relation 


(dP/dT) p 

C P = C v -T (ap/av) T 


and setting 


° v i " °Pi R 

The resulting expression for Cp is 


,( 14 ) 


c p - c Pl 


1 + ( n =± 
7 1 


e\ 2 e VP '*' 2cp 
V |”i_ e (e/T)j2 + RT 2 


,(0 /t) 


< = 


1 + 


f 2 -bp 1 cp \ 

\1— bp + 2 RT 5 / 
1, _ 2cp 

(1-bp) 2 RT 2 


(15) 


The ratio of specific heats 7 follows directly, of course, from 
equations (13) and (15) and may be written 


7 = 7 1 


1 + 


71 "I ( 

?i 


.(e/T) 


l-e 


(e/T)' 


2 + RT 2 


1 + 


2-bp 

1-bp 



_ 2cp 


( 1— bp ) 2 RT 2 * J 


1 + (rj-i) 


(D 


2 je/ t ) 

6 + 2C£ 

2 RT 2 


rw 


L 


e/m 


(16) 


A matter of special interest in the study of channel flows is the 
variation of Mach number with cross-sectional area of the channel. 

The relation between these two quantities is most conveniently deter- 
mined from the ratio of cross-sectional area at the sonic velocity 
station (the station of minimum area for any gas, as pointed out by 
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Tslen) to cross-sectional area at a station of arbitrary Mach number. 
The continuity equation yields this ratio in the form 

(*) 

With the aid of this equation and equations (l), (8), (11 ), and (12), 
the area ratio corresponding to a particular Mach number may be 
determined from a knowledge of the ideal gas specific heats, character- 
istic constants of the gas (i.e., b, c, and 0), and, for example, the 
stagnation conditions. In general, however, this computation cannot 
be carried out entirely analytically. For example, it is seen from 
equation (8) that none of the variables can be explicitly determined ; 
hence, a final solution for any single variable must be obtained 
graphically or by equivalent means. 

In order to relate the fluid properties on the two Bides of a 
normal shock wave, it is convenient to integrate equation (9)* the 
energy equation, across the wave. The resulting expression is 



Continuity of flow and conservation of momentum must also be satisfied 
across the wave. These requirements may be expressed analytically in 
the familiar forms, respectively, 

P X^1 =P 2^2 (19) 

and 

PiV-9 2 V=V*i (20) 

Equations (l8) through (20) provide the additional information 
necessary to determine the flow through a normal shock wave- Here 
again, however, final solutions for certain of the unknown flow 
parameters are most conveniently obtained graphically. 


Approximate Flows of an Imperfect Diatomic Gas 

This part of the study is concerned with first-order deviations 
of one— dimensional gas flows from the behavior of an ideal gas. With 
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this restriction, explicit solutions can he obtained, for all the 
Important flow parameters. Equation (2) will, of course, he employed 
as the equation of state. A simplification is also allowed in the 
temperature function for the vibrational contribution to the 
specific heats, since the specific heats are restricted to vary only 
in the first order from their ideal gas values. In this case only 
large values of ©/ T are considered and the Planck term may be 
approximated by 


(e/T) 


T ' [l-a 


(e/T) 


<£>' 


2 -(e/T) 

o 


Equation (8) may now be reduced to the form 

l 


In 


[G9 (?)]•■* fe-O-S 

'± + l) e-( 0 / T o) 

•To / 


. P/Po 

L(t/t q ) 2 


- 1 


-(; + 1 ) s_(e/T) (21) 


This equation can be readily solved for p/p Q in terms of T/T 0 , 
retaining only terms of the proper order, thus yielding 

i i 

- + l) e-(3/To) + [ (£) + l ] „-<e/T 0 )(T 0 /T) j (22) 

By solving for T/T 0 in a similar manner, there is obtained 


_t . .^Y 1 " 1 
To Vp 0 / 




3-27 i 


J ETo 2 L V p o 7 

1 I e -(0/T o )(Po/p) 7l_1 


- 1 




( 23 ) 
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A relation between the pressure and density ratios for isentropic 
expansion is obtained by combining equation (2) with equation (23). 
The resulting expression is 



(2*0 


Equations (22) through (24) are useful for determining the nature and 
extent of small departures from ideal gas behavior of a diatomic gas 
undergoing isentropic changes in state. 


Isentropic relations will now be obtained between pressure , 
temperature and density ratios, and the stream Mach number. In order 
to do this, it is convenient to reduce equation (12) to a form 
consistent with the first-order approximations. Performing this 
operation yields 


\ T / L ° \PoJ RT 0 2 Vpo ) 


3-S7i 


( 71-1 

■ 7 ± 


- fCvTi/ 

!_1V 

p 0 / 

J l» V 

w 

PoRTo ' 


(-&) 


_ £££o (l - Si) + ± r e -(0/T o )_ e -(0/T o )(T o /T)’l V (25) 

RT 0 ? \ T/T 0 / T 0 L J J 


where 


n = (7i 2 -37i+3) 


An expression for the temperature ratio in terms of Mach number, 
stagnation temperature and density, and gas constants is obtained by 
substituting equations (2), (22), and (24) into equation (25) and 
solving; hence 
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where 




It may be noted that equations ( 2 6) , (28) , and (30) can be rearranged to give the Mach number in 
terms of T/T 0 , p/p 0 , and P/P 0 , respectively. 
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Hence, substituting equations (26), (28), and (35) into equation (17) in the above form yields 
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Substituting equations (2), (19) , and (20) into equation (38), dividing through by equation (35) 



and solving for the -temperature ratio yields the expression: 
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In order to obtain the relation for the density ratio Oz/Qx, it is convenient to write 
equation ( 39 ) i n the form 
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is the sum of all the first-order correction terms. The pressure ratio across a normal 
is now obtained, by combining equations ( 19 ), (20), (35), and (44) to provide the 
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upstream of the wave with the three identities 
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Po 2 _ ( P 2 Y Piftoi > 
P 0l \ ?! /\Pi 7 *^/ 


!z* Y ^ Y Pl/Po O 

p Ql \ P x Ap 2 /p 02 / 

and 

To^ _/ N 

T 0l \ T X At 2 /To 2 / 


(^9) 


(50) 


(51) 


( in which all terms of the right-hand members are given in previously 
developed expressions. 


It is important to note that the normal shock expressions are 
applicable to the study of oblique shock flows. For example, 
equations ( 39 ) through (46) may be employed directly by substituting 
M x sin 3 for M* wherever It occurs. In order to determine Ms, 
however, it is convenient to have a relation between the wave angle 3 
and the stream deflection angle 6. This relation may be obtained by 
combining the momentum, energy, and continuity equations for flow 
through an oblique shock (reference 9 ) to yield the expression 


ff - 1 


_ PjJi 


- ( 
i v 


,cot 3 cot 6 


tt) 


which may be written 


cot 5 = tan 3 


L (1-n) lPi7FTPI 


- 1 


(52) 


M 2 can then be calculated from the expression 


Ma = (t £\ (1 - g P + % ) 

Mi VTa/i Vp 2 A sin (3-6) 


(53) 


which is obtained from the continuity equation. In this equation 
the ideal gas values for the temperature and density ratios are 
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those corresponding to the initial Mach number Mi and the wave 
angle p. Knowing Mb from equation (53), the remaining normal 
shock expressions, equations (49) through (51), can also be applied 
to oblique shock flows. 

It is observed that, in general, all expressions obtained in 
the preceding approximate analysis differ from the corresponding 
ideal gas equations only by small correction factors containing the 
Mach number, temperature and density, and constants of the gas. As 
pointed out in the introduction, this simplifies the investigation 
of one-dimensional imperfect gas flows, provided effects of the 
imperfections are not large. 


DISCUSSION AND CONCLUSIONS 

The complex manner in which gaseous imperfections influence 
general one-dimensional-flow processes precludes the detailed consid- 
eration, here, of more than two flows of aerodynamic interest. 

Before discussing these flows, it may be worthwhile to establish more 
definitely the conditions under which imperfections in air can alter 
its dynamical behavior from that predicted on the assumption that 
ideal gas laws are obeyed. These conditions are easily deduced from 
the equations employed in the analysis and a knowledge of the char- 
acteristic gas constants . 3 They may be summarized as follows : 

1. The occurrence of temperatures in excess of 800° R, in which 
case the specific heats change appreciably. 

2. The occurrence of densities in excess of 10 times sea level 
atmospheric, in which case the molecules occupy an appreciable 
fraction of the available volume. 

3. The occurrence of sufficiently high, pressures and low 
temperatures to cause the intermole cular forces. to be appreciable. 

One additional observation should be made concerning the inte- 
grated effect of these phenomena on a particular flow process. This 
is that the several imperfections frequently counteract each other 
(note, e.g., the state equation and equations (24), (26), and (27)), 

3 The values chosen for the characteristic gas constants, and the 
method by which they were determined, are included in Appendix B. 
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and. thus alter certain flow parameters only slightly from their ideal 
gas values . 

The special cases of flow through normal and oblique shocks in 
free air (NACA standard atmospheric conditions at sea level) will 
now be considered. The pressure and temperature ratios across a 
normal shock have been calculated for Mach numbers from 1 to 10 
using equations (l), (ll), (18), (19), and (20). In this case the 
terms containing b and c are extremely small (they could, in 
fact, be neglected), and it is sufficiently accurate to use the 
values for these constants given in Appendix B. 

The pressure rise so determined is shown in figure 1 as a funo- 
tion of Mach number, and it is seen that it differs by less than 
6 percent from the ideal gas value up to the maximum Mach number of 
10. This behavior might be anticipated from the ideal gas equation 
for the pressure ratio (Appendix A) which is relatively insensitive 
to variations in 7 at high Mach numbers . The corresponding tempei>- 
ature rise across a normal shock is illustrated in figure 2. Here 
it is observed that the effect of increase in specific heats appreci- 
ably reduces the temperature ratio below the ideal gas values at the 
higher Mach numbers 5 . It may be deduced from figures 1 and 2, that 
the density rise will be considerably above the ideal gas value at 
the higher Mach numbers. These results must be qualified at Mach 
numbers in excess of 7, for under these circ ums tances the temper- 
atures downstream of the wave appreciably exceed 5000° R which is 
sufficiently high to cause molecular dissociation and thus possibly 
to alter the flow. 

In view of the fact that the changes with temperature of the 
specific heats of air do not strongly influence the pressure rise 
across a normal shock, up to Mach numbers of 10, it follows that 
corresponding pressure effects in flow through oblique shocks will 
also be small, provided the Mach number of the component of velocity 
normal to the wave is below 10. This suggests, within the limita- 
tions observed, that the effect of variable specific heats on those 
aerodynamic characteristics essentially independent of viscosity may 
be small for some shapes. (Epstein* s findings in reference 11, 
concerning hypersonic flows about bodies, pertain to much higher 
Mach numbers that those considered here.) The extent to which this 

4 

All ideal gas data are obtained from the tables of reference 5. 

5 This result is consistent with Wood's findings concerning surface 
temperatures at high Mach numbers. (See reference 10.) 
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is correct has been investigated for the simplest shape, a flat 
plate. Pressures on the upper surface of the plate were determined, 
assuming Prandtl-Meyer flow about the leading edge, and pressures 
on the lower surface were found with equations (45) (replacing M x 
with Mx sin 3) and (52). These data and the corresponding lift and 
pressure drag coefficients of the plate were calculated for Mach 
numbers of 10 and 20, and for angles of attack from 0° to 24° and 
13°, respectively. A comparison of the coefficients with those 
obtained assuming ideal gas flow over both the lower and upper 
surfaces showed that the effect of variable specific heats is very 
small. (E.g., the force coefficients are smaller in the case of the 
imperfect gas but never by more than 3 percent, even at the highest 
angles of attack.) This result is not justification for assuming 
ideal gas behavior to calculate the high Mach number aerodynamic 
characteristics of arbitrary shapes, other effects such as specific 
heat lag being neglected; but it does indicate that the assumption 
is justifiable for certain shapes on which the important pressure 
forces are obtained through oblique or normal shock compression. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 
IDEAL GAS EQUATIONS 

The following are ideal gas equations used in the analysis of 
approximate flows of an imperfect diatomic gas: 


Isentropic Flow Equations 



Normal Shock Equations 


„ 2 (7i-l) M! 2 + 2 

Ms r- - 

1 27iMi 2 - (/i— 1) 

‘s' \ _ 27jMi g - (71-1) 

VPxA 7i+l 
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APPENDIX B 

DETERMINATION OF CHARACTERISTIC GAS CONSTANTS FOR AIR 

Numerical values for the characteristic gas constants b «r>H c 
may be determined from a knowledge of the critical pressure and tem- 
perature of a gas. Values determined in this mann er are, however, 
generally .less accurate in a particular range of temperatures A nd 
pressures than is desired. Consequently, it is often advantageous to 
fix b and c such that theoretical and experimental values of the 
primary parameters o p , 7 and compressibility (i.e., p/p) are in 
good agreement in the range of interest. This method is the more 
suitable of the two for most aerodynamic studies as the variations 
in pressure and temperature are limited; hence, it has been anployed 
in this paper to determine b and c, as well as 9, for air. 

It is sufficient for the illustrative applications presented here 
to determine the former constants only for the approximate flow 
equations. Separate values of 9 are found for the exact and approxi- 
mate Planck terms. 

As pointed out previously (note the development of approximate 
expressions for c p , c v > and 7) , it is a simple matter to obtain the 
characteristic gas constants appearing in the approximate flow equa- 
tions. In this example, c and 9 were chosen first to give 
agreement between experimental and theoretical values of c p and 7, 
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the latter values "being calculated from equations (33) and (3*0, 
respectively. Data on the variation of these quantities with 
pressure and temperature at high temperatures were obtained from 
references 12 through 19. 0 A comparison of these data with theory, 
setting c = 2.25 x 10® °R,ft 5 /slug,sec 2 and 9 = 58OO 0 R, is shown 
in figures 3 end 4 for pressures of 0 and 144 X 10 pounds per 
square foot absolute. The agreement is observed to be good up to 
temperatures of 3000° R. Using equations (15) and (l 6 ) to calculate 
Cp and 7 , respectively, and a value of 9 - 5500° R, it is seen that 
at zero pressure excellent agreement with the correlated data is 
obtained up to 5000° R. The approximate theory was then checked at 
low temperatures with experimental values of 7 given in reference 
20. At temperatures above liquefaction and pressures up to 25 
atmospheres, the difference between theoretical and experimental 
values of 7 was found not to exceed 3 percent. 

The molecular size constant b was chosen to yield good agree- 
ment between theoretical and experimental values of compressibility 
of air. The theoretical values were calculated from equation (2), 
and the experimental data were obtained from reference 21. In 
general, it was found that using b = 0.485 ft /slug and the previ- 
ously determined value of c equation (2) predicts the compressi- 
bility of air to within 2-percent error at pressures from 0 to 150 
atmospheres, and temperatures from liquefaction to 850 0 R. (The 
deviation from ideal gas behavior is small above this temperature 
at the pressures indicated. ) 
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Figure f. — Variation of pressure ratio across a normal shock 

wave with Mach number. (Initial state- NACA standard 
atmosphere at sea level.) 
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Figure Variation of temperature ratio across a normal 
shock wove with Mach number. (Initial stote-NACA standard 
atmosphere at sea level.) 
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Figure 3 .— Variation ofC p with temperature at absolute pressures of 0 and 144 x 10 
lb / sq ft ebs. 
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Figure 4.- Variation of 7 with temperature at absolute pressures of O and 144 x/o 3 
ib / sg ft obs. 




